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ABSTRACT 

We have conducted an H I 21 cm emission- line survey using the Parkes 20cm multibeam instrument 
and the Australia Telescope Compact Array (ATCA) of six loose groups of galaxies chosen to be 
analogs to the Local Group. The goal of this survey is to make a census of the H I-rich galaxies 
and high-velocity clouds (HVCs) within these groups and compare these populations with those in 
the Local Group. The Parkes observations covered the entire volume of each group with a rms Mhi 
sensitivity of 4-lOxlO 5 M0 per 3.3 kins' 1 channel. All potential sources detected in the Parkes 
data were confirmed with ATCA observations at ~2' resolution and the same Mjji sensitivity. All the 
confirmed sources have associated stellar counterparts; no starless H I clouds-HVC analogs-were found 
in the six groups. In this paper, we present a description of the survey parameters, its sensitivity and 
completeness. Using the population of compact HVCs (CHVCs) around the Milky Way as a template 
coupled with the detailed knowledge of our survey parameters, we infer that our non-detection of 
CHVC analogs implies that, if similar populations exist in the six groups studied, the CHVCs must be 
clustered within 90 kpc of group galaxies, with average M H i< 4xlO 5 M at the 95% confidence level. 
The corollary is that the same must apply to Milky Way CHVCs. This is consistent with our previous 
results from a smaller sample of groups, and in accordance with recent observational and theoretical 
constraints from other authors. These results confirm that there is very little neutral matter around 
galaxies, and that any substantial reservoir of baryons must be in other phases. 
Subject headings: galaxies: formation — intergalactic medium — Local Group 



1. INTRODUCTION 

Th e majority of galaxies reside in groups (jTullvl 
1987). Loose groups of galaxies, such as the Lo- 
cal Group in which the Milky Way is located, are 
a collection of a few large, bright galaxies and tens 
of smaller, fainter ones. The large galaxies are typ- 
ically separated by a few hundred kiloparsecs and 
spread over an extent of approximately a megaparsec. 
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Groups of galaxies are the building blocks of galaxy 
clusters, feeding gas-r i ch sp ir al galaxies into the clus- 
ter (e.g. [Burns et al.l 11994 iBravo- Alfaro et alJ 120001 : 
iGonzalez et al.1 120051 ). Loose groups themselves are 
still c ollapsing and not virialized (|Zabhidoff fc Mulchaevl 
Il998h . The gas-rich galaxies within loose groups are also 
still form ing as they accrete their dwarf g alaxy satel- 
lites (see iFreeman fe Bland-Hawthornl I2002L and refer- 
ences therein) and their reservoir of gas (e.g. iBlitz et ail 
[1999I l Some of this reservoir may be in the form 
of neutral gas and possibly related to the "high- 
velocity clouds" (HVCs) seen around the Milky Way (see 
iWakker fc van Woerdenl \l99% \ van Woerden et all |2004| . 
and references therein for a review). 
HVCs are clouds of neutral hydrogen (H I) discovered 
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approximately 40 ye ars ago (iMuller et al .11 1963ft covering 
up to 37% of the sky ([Murphy et al.l ll995: Loc kman et al.l 
2002) with velocities inconsistent with simple Galactic 
rotation and in excess of ±90 kms" 1 of the Local 
Standard of Rest (jWakker fc van Woerde"nlll997ft They 
lack associated stellar emission dSimon fc Blitd [2002; 
IWillman et"aTll2002t ISiegel et al.ll2005t ). and so we have 
no direct measure of their distances, and, therefore, their 
masses. As such, we can not easily discriminate between 
their possible origins. Nevertheless, HVCs likely repre- 
sent a variety of phenomena. 

Some HVCs may be related to a g alactic fountain 
(|Shapiro fc Field! [19761 iBregmanl Il98l . Other HVCs 
are certainly tidal in origin: the Magellanic Stream is 
the most obvious of these features, formed via the tidal 
interactions between the Milky Way, La rge Magellanic 
Clou d, and Small Magellanic Cloud (e.g. iPutman et al.l 
120021) . Other HVCs may b e related to the Sagittarius 
dwarf (|Putman et all 120041) . Some HVCs may e v en be 
satell ites unto themselves ( Lockmaril 120031 ). lOortl (|1966l 
1970) originally proposed that HVCs may be infalling 
primordial gas associated with the formati on of the Milky 
Way. Complex C may be such an example ( Wakker et al.l 
1999: iTripp et al.ll200l iGibson et al1l200lL 

While the idea of associating HVCs with galaxy for- 
mation is not ne w, recently t his hy poth esis has attracted 
more attention. IBlitz et all (|1999f ) and lBraun fc Burton! 
(1999) suggested that the HVCs and the compact HVCs 
(CHVCs), respectively, may contain dark matter and 
could be related to the small dark matter halos predicted 
to exist in large numbe rs by cold dark matter models of 
galax y formation (e.g. iKlypin et al.l Il999t iMoore et all 



1999)" As originally conc eived bv IBlitz et alJ (|1999f ) 
and lBraun fc Burtonl ([1999ft . HVCs or the subset of the 
most compact HVCs (CHVCs) have distances ~1 Mpc, 
and M^f/~1O 7 M0. Current models associ a ting C HVCs 
with dark matter halos by Ide Heij et al.l (|2002bft sug- 
gests the CHVC distribution has a Gaussian distribu- 
tion about the Milky Way and M31 with D-150-200 kpc 
and Mgj~1 5 - 5 ~ 7 M(7). A similar dist ance distribution is 
suggested by lMaller fc Bullockl (|2004l ) who propose that 
the HVCs are the products of the cooling, condensing 
hot halo. Such clouds, however, would not be associated 
with dark matter halos. 

There have been many detections of possible HVCs 
around other galaxies. In these systems, we can easily 
determine the masses and separation of the clouds. We 
can clearly see examples of extragalactic H I clouds as- 
sociated with galactic fountains, tidal debris, and even 
very faint dwarf galaxies, but none to date can be un- 
ambiguously identified as being associated with galaxy 
formation or dark ma tter halos. The H I Rogues Gallery 
ijHibbard et alJl200ll) 3 contains many examples of such 
H I clouds. 

ISchulman et alJ (|1994l ) identified high-velocity wings 
on integrated H I profiles from a sample of face-on galax- 
ies. The presence of these wings was well-correlated with 
an increased star formation ra t e in t he galaxy, and so it 
is likely that ISchulman et al.l (I1994D identified H I as- 
sociat ed wi th a galactic fountain iKamphuis fc Sancisl 
(|1993ft and IBoomsma eFaTI (|2004ft have identified high- 
velocity H I associated with NGC 6946. Much of this 

3 available online at http://www.nrao.edu/astrores/HIrogues/ 



gas is coincident with H I holes and is predominately 
coincident with the optical disk. As NGC 6946 is pro- 
lifically forming stars, and has had eight supernovae 
in the past century, it seems likely that most of this 
gas must arise from a galactic fountain. O ther galax- 
ies, such as NGC 891 (iSwaters et all fl99l . NG C 253 
(IBoomsma et al.ll2005ft. N GC 4559 dB arbieri et al.ll2005ft . 
NGC 2403 (|Fraternali et al.ll2002D . also show signatures 
of extra-planar H I possibly associated, at least in part, 
with a galactic fountain. 

We know of many examples of gaseous tidal debris 
without optical counterparts around individual galax- 
ies. The Leo Ring is a 200 kpc diameter ring of H I 
orbiting within the M 96 Group ([Schneider et al.|[l983l : 
lSchneided[l985ft . While it may have a tidal origin, as i t 
appears to be interacting with M 96 (|Schneidedll985ft . 
its orbital period of 4 Gyr and its coherence over such 
a large area is a challenge for such an explanation. In- 
stead, it may be that this ring is a primordial remnant 
of the formation of the group rtSchneiderl 119851 ). A -50 
kpc diameter ring (—100 kpc in length), a similar fea- 
ture to the Leo Ring b ut on a much sma ller scale, is 
seen around NGC 4449 (|Hunter et al.lll998ft . While this 
could be explained as o ngoing accretion of primordial 
gas ([Hunter et alJ Il998h . the structure of this ring (if 
not its presence) can be explain ed via a recent inter - 
action with a nearby companion (|Theis fc Kohkj f2001). 
IC 10 has an extended H I distribution that is counter- 
rotating with respect to the main body of the galaxy, 
plus a long streamer e xt ending away from the g alaxy 
(|Wilcots fc Milled 119981 ). IWilcots fc Milled (|1998f i sug- 
gest this H I represents ongoing infall of primordial mate- 
rial, and the absence of a nearby perturber, despite resid- 
ing in the Local Group, provides circumstantial support 
for this scenario. NGC 925 represents our final exam- 
ple of intergalactic H I which may be primordial or tidal 
in origin. NGC 925 has a -10 7 M^ H I cloud attache d 
to the main galaxy by a streamer (|Pisano et al.lll998ft . 
NGC 925 does reside in a group, but the nearest galaxy 
to NGC 925 is 200 kpc away. Is this H I cloud then pri- 
mordial gas, the debris from an old tidal encounter, or a 
d isrupted dwar f galax y? 

iPisano et alJ (|2002ft searched for H I clouds around 41 
quiescent, isolated galaxies and found only gas-rich dwarf 
galaxies. The absence of starless H I clouds around iso- 
lated galaxies suggests that the presence of such clouds 
in denser environments (such as those described ab ove) 
may be the result of interactions (jPisano et~aT1l2002T ). To 
date, all claims of intergalactic H I clouds without asso- 
ciated stars, or "dark gal axies" (e.g. H I 1225+0146 
iGiovanelli fc Havnesl |1989| ) have turned out to be, on 
closer inspection, eithe r low surface brightness galaxies 
(|McMahon et al. | ll990ft or tidal de bris connected with a 
bright galaxy (|Bekki et all l2005alTbl ). These detections 
tend to be in dense environments, where tidal interac- 
tions are more probable. As such, if we want to find 
HVC analogs associated with galaxy formation, it is im- 
portant to search for such H I clouds, using very sensitive 
instruments, in low density environments similar to the 
Local Group. 

There have been many previous surv eys of galaxy 
group s searching for starless H I clo uds. lLo fc Sargent 
mm, [K raan-Korteweg et al.l |l999), and iKovac et alJ 
(2005) surveyed the entirety of the Canes Venatici I 
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group. They found many new d warf g alaxies, but no 
HV C analogs. IZwaan fc Briggsl (|2000h . IZwaanl ([20011 ) 
and Ide Blok et alj (|20o3 T conducted a sparse survey of 
nearby groups, again finding no H VC analogs. The 
HIPASS survey (|Barnes et al.l [200lh covered the entire 
south ern sky, and, whi le some H I clouds were detected 
(e.g. iRvder et al.ll200ll ), these are likely tidal in origin. 
The deepest survey of a group to date was the HIDEEP 
survey of the Cen A group, with a baM^i limit of 
2 x 10 6 M Q , but all of the new H I detections in th is sur- 
vey had optical counterparts (jMinchin et alj |2003). The 
Cen A group is dominated by a large elliptical galaxy, 
so it may not be a good analog to the Local Group, un- 
like those studied by the other collaborations, and, hence, 
perhaps not the ideal location to sea rch for HVC a n alogs . 

Our survey, as first reported in iPisano et alj (2004, 
hereafter Paper I), seeks to provide the deepest obser- 
vations of the entirety of six loose groups, all analogous 
to the Local Group, to find, or place strong limits, on 
the presence of CHVC analogs around these galaxies. 
Paper I reported the initial results of this work for the 
first three groups. In this paper (Paper II), we will dis- 
cuss the survey properties in greater detail, including its 
sensitivity and completeness, and present the results of 
our search for CHVC analogs for all six groups We use 
our mea sure of completeness to improve our model from 
iPaper II and use it to determine the implications of our 
non-detections for the distribution of CHVCs in these 
groups and around the Milky W ay. In the final paper of 
this series (jPisano et al.l [20061 hereafter Paper III), we 
will examine the properties of the galaxies in these loose 
groups and how their properties compare to the popula- 
tion of the Local Group. 

This paper is structured as follows. In Section [2] we 
discuss our sample selection and the basic properties of 
the groups we observed. The Parkes observations and 
reductions are presented in Section [3j followed by a dis- 
cussion of source-finding in Section!?] The confirming in- 
terfcrometric observations and reductions are discussed 
in Section [5] In order to determine the implications of 
our failure to detect any HVC analogs, we must have a 
good understanding of the completeness of this survey; 
this is presented in Section [6] The results of our survey 
are briefly summarized in Section [7] The implications for 
the distribution of CHVC analogs and the nature of the 
HVCs around the Milky Way are discussed in Section 
and we conclude in Section [9] 

2. SAMPLE SELECTION 

We selected nearby («gsr < 1000 kms" 1 ) loose groups 
of galaxies that are analogs to the Local Group, which 
contain only spiral and irregular galaxies. For these 
nearby groups, we can obtain high spatial resolution and 
can detect very low-mass H I clouds. Groups that may 
be confused with Galactic or Local Group H I emis- 
sion (wgsr < 300 kms -1 ) were avo ided. In the "Lyon 
Groups of Galaxies" catalog CLGG: [Garcial fl993h eight 
groups in the Parkes declination range (S < 0°) matched 
these criteria. We chose to observe five groups of galax- 
ies from this catalog: LGG 93, LGG 106, LGG 180, 
LGG 293 and LGG 478. In addition, we chose a sixth 
gr oup with sim ilar properties from the group catalog 
of IStevensl (120051) . which is solely composed of HICAT 
galaxies (jMever et alj 12004) that we call "the HIPASS 



Group" . Properties of these groups are listed in Table [T] 
and dis cussed below. We assume an H = 72 kms -1 
Mpc -1 ([Spergel et al.ll2003t ), and calculated the distances 
to the groups after correcting their heliocentric velocities 
for a multi-attractor velocity flow model of the local uni- 
verse (Masters 2005, Masters et al. 2007, in preparation) . 
Note that as a result of this new m ethod, t he distances 
used here are different from those in lPaper^ . The uncer- 
tainties on these distances are ~2 Mpc (Masters, 2006, 
priv. comm.). 

3. PARKES OBSERVATIONS AND REDUCTIONS 
3.1. Multibeam Observations 

The six groups were observed, in sessions of 
approximately 10 nights each, between October 
2001 and June 2003 using t he 2 0cm multibeam 
instrument (|Stavelev-Smith et al.l Il996h on the Parkes 
64m radiotelescope. The multibeam has 13 hexagonally- 
packed receivers. Each individual half-power beamwidth 
is 14.3' in diameter and is separated by two beamwidths 
from its neighbor. The whole system spans an on-sky 
diameter of 1.7°. For LGG 93 and LGG 180, we ob- 
served with an 8 MHz bandwidth using only the inner 
seven beams. The velocity resolution for these observa- 
tions was 1.65 kms -1 . By August 2002, new 16 MHz 
filters were available for all 13 beams, and were used for 
the observations of the remaining four groups. The 16 
MHz bandwidth observations had a velocity resolution of 
3.3 kms -1 . Only data taken at night was used to avoid 
solar interference which can both add noise to the data 
cubes and result in spurious detections. 

All groups were observed by scanning the multibeam 
across the group along great circles in both the right as- 
cension and declination directions at a rate of 0.5 degrees 
per minute. With data being recorded every 5 seconds, 
each beam area is sampled almost six times in the scan 
direction. For observations with the inner seven beams, 
we rotated the multibeam by 19.1° such that the half 
power points of each beam touch. For observations using 
all 13 beams, the rotation angle was 15° to approximate 
the same coverage. Each scan was offset by 13' so that 
the field is uniformly sampled perpendicular to the scan 
direction by each beam. We cycled through all the scans 
covering the full area of the group, in both directions, 
until the desired sensitivity was reached. The charac- 
teristics of our observations for each group are listed in 
Tabled 

3.2. Calibration 

A default calibration based on past measurements is 
applied to the data as it is written to disk. This default 
calibration was determined for the HIPASS setup with 
64 MHz bandpass and a central frequency of 1384 MHz. 
We regularly (usually once a day) observed Hydra A to 
check the calibration using our observing configuration, 
and used our derived average calibration factors instead 
of the defaults in our reduction process. The derived 
factors were typically within 10% of the default values. 

3.3. Reductions 

All data were reduced and gridded using the livedata 
and GRIDZILLA packages in AIPS++. The details of how 
these packages work and were used for reducing HIPASS 
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data are discussed in iBarnes et al.l (|2001[ ). We used a 
two-step iterative process in reducing our data. 

livedata was used to remove the bandpass both tem- 
porally and spectrally. The bandpass was fit in the scan 
direction with a first order polynomial, where outlying 
points (from the source or interference) were clipped at 
the 1.5 a level over the successive iterations. This is 
equivalent to subtracting and normalizing by an 'off' 
position for a standard 'on-off' observation. A first or- 
der spectral baseline was then fit to the resulting spec- 
trum at each position of each individual scan. Like for 
the HIPASS reductions, a 25% Tukey smoothing func- 
tion was applied to damp out the ringing resulting from 
strong line sources (like the Galaxy). The reduced scans 
were then combined into a single data cube using the 
program GRIDZILLA. GRIDZILLA determines which spec- 
tra contribute to which pixel and with what weight they 
contribute, while discarding outlying spectra. We used a 
mean weighting scheme while clipping the brightest 2% 
of the spectra. Since there were a few hundred spec- 
tra contributing to a given pixel, this resulted in a very 
small loss in signal-to-noise, while removing most of the 
effects of interference. We were able to use mean grid- 
ing as compared to the median griding used for HIPASS 
because our data were all taken at night eliminating con- 
tamination from solar interference, but some other forms 
of interference are visible in our data cubes. This also 
results in an improvement in the noise. 

At this point, the cubes were searched for sources 
(more on this in Section |4|). All sources were then 
masked, and the reduction process was repeated with 
a few small differences. In livedata, a second order 
polynomial was fit in the time (spatial) direction with 
2a clipping over three iterations. A second order polyno- 
mial was also fit and removed from the spectral baselines. 
Because the sources are masked, these higher order fits 
yield a flatter residual baseline while lacking large nega- 
tive sidelobes around bright sources that can hide nearby 
weak sources. Furthermore, this process yields more ro- 
bust measures of the H I spectrum of all sources. The 
resulting noise level and mass sensitivities are listed in 
Table H 

4. SOURCE FINDING 

Once the data cubes for each group were complete, 
three teams of authors (DJP, DGB, and BKG & VAK in 
tandem) searched each cube by eye for detections. Any 
source identified in two or more searches was considered 
"real" and was a candidate for follow-up confirmation 
with the Australia Telescope Compact Array 4 . Each 
cube had simulated sources inserted into the data by a 
third party (Dr. M Zwaan) to provide a measure of the 
completeness of our survey as a function of linewidth and 
integrated flux. This process was done i n the same fash- 
ion as for HIPASS (| Zwaan et al.l I2004D . As described 
above, after all the possible detections were identified, 
the final, masked cubes were made and re-searched to 
identify any new sources-none were found. 

5. ATCA OBSERVATIONS AND REDUCTIONS 



We used the ATCA to observe all 105 of the 112 Parkes 
detections in and behind the six groups to confirm the 
reality of all our detections, to uniquely identify optical 
counterparts, and resolve any confused sources. The re- 
maining seven galaxies were behind one of the groups 
and previously detected in HIPASS. Of the remaining 
105 possible sources, 15 had previously been observed 
with the ATCA or VLA. We observed the remaining 90 
between October 2002 and March 2005 using a compact, 
750m configuration for observations of all of the groups 
to obtain a nearly circular beam of ~l-2'. Because of the 
near equatorial position of LGG 293, it was the excep- 
tion; we observed its galaxies with the H214C configura- 
tion (a configuration with antennas on the ATCA's North 
Spur with a maximum baseline of 214m) for resulting res- 
olution of -2-3'. We used a bandwidth of 8 MHz (~1700 
kms -1 ) with 512 channels of 3.3 kms -1 after Hanning 
smoothing. We observed 1934-638 as a primary flux and 
bandpass calibrator at the beginning or end of each day's 
observing run, and made hourly observations of fainter, 
closer radio sources for phase calibration. To confirm 
all our detections, every source was observed with the 
ATCA to at least the same sensitivity as in the Parkes 
data, typically about 4 mJy/beam. The majority of the 
observations occurred at night to minimize the effects of 
solar interference. 

The calibration and reduction of the data was done in 
the standard way using MIPJAD. Data were flagged by 
hand to remove interference, both terrestrial and solar. 
Generally very few data needed to be flagged. The data 
for each day were separately reduced and calibrated. The 
line-free region in each cube was identified and the data 
were continuum subtracted. The line data from each 
night were then combined to make a data cube using a 
robust (a robustness of two) or natural weighting scheme. 
The cubes were then CLEANed down to 2 — 3 times the 
RMS noise level, whichever was sufficient to remove all 
CLEAN artifacts. Moment maps of the total intensity 
and velocity field (moments and 1) were made in AIPS 
with 2 — 3cr blanking over the velocity range where emis- 
sion was observed. 

As our groups are nearby and composed of large, bright 
galaxies, many of the galaxies in our groups had been 
previously observed with the VLA and the ATCA. We 
have taken these data from the relevant on-line archives 
and re-reduced the data in the same manner described 
above. This archival data generally had a better sen- 
sitivity than our Parkes data, with velocity and spatial 
resolution comparable to our ATCA 750m array obser- 
vations. 

All cubes were inspected by eye for H I emission from 
the Parkes detection and for additional emission from 
HVC analogs or gas- rich dwarf galaxies. Of the origi- 
nal 112 Parkes detections, only five were not confirmed 
with the ATCA or VLA. An additional four dwarf galax- 
ies, behind the target groups, were resolved in the inter- 
ferometric data. No HVC analogs were detected in the 
ATCA or VLA data, despite the higher spatial resolu- 
tion and equivalent mass sensitivity as compared to the 
Parkes d ata. All o f these data will be presented in more 
detail in iPaper Till . 



4 The Australia Telescope Compact Array (ATCA) is part of 
the Australia Telescope, which is funded by the Commonwealth of 
Australia for operation as a National Facility managed by CSIRO. 
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In order to obtain an accurate census of the H I-rich 
objects in the six loose groups, it is essential to have a 
measure of the completeness of our survey, particularly 
the completeness as a function of linewidth, W20 and in- 
tegrated flux, Si n t, as we derive mos t propert ies of the 
galaxies from these two quantities. In lPaper I . we made 
the blanket assumption that our survey was 100% com- 
plete for sources with an integrated flux brighter than 
lOer, and 100% incomplete for fainter sources. In this 
paper, we refine our measures of completeness. 

The completeness was evaluated in three ways. First, 
we examine the distribution of the number of sources 
(real and simulated) as a function of W20 and Si n t inde- 
pendently. Second, we examine the distribution of simu- 
lated sources as a joint function of W20 and Si n t ■ Finally, 
we compare our data with the completeness function for 
HIPASS. All methods have their advantages and draw- 
backs. 

Figure [1] shows the number of real and simulated 
sources (both detected and undetected) as a function of 
W20 and Si n t ■ Because this plot represents the composite 
completeness for six fields with slightly different sensitiv- 
ities and channel sizes, we must find a way to scale the 
parameters to a common scale. We use the following 
formulae to do so: 



SNR V 



W20 
AV 

Sint 



(1) 

(2) 



where AV represents the channel size in kms -1 and a 
represents the noise in a single channel. A 2 o and SNR int 
are in units of channels and integrated signal-to-noise ra- 
tio as compared to the theoretical limit if the noise is 
Gaussian. The completeness is just the ratio of the num- 
ber of detected simulated sources to the total number, 
and this has been applied to the data in the left panels 
(as shown by the open circles) . It is evident that the drop 
in the number of real sources as a function of SNRint is 
purely due to incompleteness and does not reflect the 
true flux distribution of our sources. 

For A20 the completeness is much more uniform, as 
seen from the simulated sources. In fact, the incomplete- 
ness appears to have little dependence on the velocity 
width of the source. The drop in the number of real 
sources as a function of velocity width is, again, not just 
due to incompleteness. We could have smoothed the data 
and searched the cubes over a range of velocity resolu- 
tions. In principle, this would make it easier to detect 
sources whose velocity widths best match our resolution 
elements. Since we did all our searching by eye, we chose 
not to do this due to the labor involved. 

Figure [2] shows the properties of all simulated sources, 
and those which were identified by more than two au- 
thors. As there were only 10 simulated sources per 
cube, we have used Equations [T] and [2] to place all the 
sources on a single plot. It is evident from Figure [2] that 
sources fainter than 5a are not detected regardless of 
their linewidth. Above this flux level, the presence of 
ripples in the spectrum, caused by standing waves from 
broadband interference can degrade our ability to detect 
broad sources. Even for sources brighter than 10er we are 
only 76% complete. This is due to broad line sources with 



low peak fluxes, only two of ten sources with linewidths 
smaller than ten channels and fluxes above 10tr are un- 
detected, but is also the result of small number statis- 
tics. Note also that we have no simulated sources with 
linewidths smaller than three channels. Such sources 
would be indistinguishable from narrow band radio inter- 
ference, and, as such, highly unreliable. While we could 
use this figure directly to correct our survey for incom- 
pleteness, there are relatively few simulated sources here 
(only 60) which makes it difficult to accurately assess the 
completeness throughout this parameter space. As such, 
we will compare this to the completeness measured for 
H IPASS. 

IZwaan et al.l ()2004f ) used large numbers of simulated 
sources inserted into the HIPASS cubes to characterize 
the completeness of that survey and parameterize it as a 
function of SVnt, S pea k, and W2o- Again, we scale both 
Si n t and W20 by the HIPASS noise and channel size; 
the result is shown in Figure [3] The simulated sources 
inserted into our data are also shown to illustrate that 
these two measures of completeness are roughly consis- 
tent. Note that the completeness is not proportional to 
VW, but degrades more rapidly than that. For our anal- 
ysis, we will use the HIPASS completeness function (as 
scaled for each group), but assume 100% incompleteness 
for sources with SNRint < 5o\ While this appears to 
be reasonable, the reader should remember that HIPASS 
had different observing parameters (e.g. channel size and 
sensitivity) and observing technique (pure Declination 
scans vs. basket- weave) and a different source-finding 
approach (automated vs. manual). This completeness 
function may, therefore, not be entirely accurate for our 
survey. Nevertheless, it is probably a reasonable approx- 
imation and what we will use for the remainder of this 
paper. 

For reference, our 5a Mhi detection limit for a source 
with a velocity width of 30 kms -1 ranges from 0.5- 
2x 10 7 M© (depending on the group). The lowest mass 
detection in any of our groups has Mhi= 1 X 10 7 M Q . At 
these level, our survey is quite incomplete. 

7. SURVEY RESULTS 

A total of 111 galaxies were detected by the combined 
Parkes/ATCA survey. Of these, 63 are in the groups 
while the remaining 48 are backgro und galaxies. All of 
the galaxies in the optically-selected [Garcia (|199 3l) group 
catalog were det ected in H I as we re all galaxies in the 
HICAT catalog (jMever et al.ll2004D . These 63 galaxies 
represents a doubling of the number of galaxies in the 
optically-selected groups and a 60% increase over the 
number of HICAT galaxies in these groups. 

All the detected galaxies have optical counterparts, 
either cataloged in NED or visible in the Digital Sky 
Survey. There were no intergalactic H I clouds without 
optical counterparts detected-no HVC analogs detected. 
The properties of the newly detected group galaxies are 
consistent with them being late- type spiral, irregular, 
and dwarf irregular galaxies; only one background galaxy 
is classified as an early-type and is relatively gas-poor. 
The properties of all our detections and the ensemble 
properties of the grou ps will be presented and discussed 
in detail in lPaper III] . 

8. DISCUSSION 



G 



iPaper J contained the limits on the distances to the 
compact high- velocity clouds (CHVCs) based on our non- 
detection of any analogs in three of the loose groups 
we surveyed and a simple model for the distribution of 
CHVCs. In this paper we will use our improved char- 
acterization of the complet eness of our survey, in con- 
cert with our model from IPaper J . to place limits on 
the CHVC distribution for our entire sample of six loose 
groups. Furthermore, we will compare the limits for both 
a Gaus sian distan ce distribution and a iNavarro et al.l 
(NFW; [l996Lfl99l distribution of CHVCs, and examine 
the robustness of our models in light of their potential 
weaknesses. 

8.1. A model for CHVCs 

Our simple model was originally described in iPaperll . 
That description is presented again here. This model ba- 
sically assumes that the distribution of CHVCs around 
the Milky Way is representative of the distribution 
around other galaxies. As such the non-detection of 
CHVCs within our sample of six groups provides a con- 
straint on their distribution around the Milky Way. 
There a re 270 CHVCs in th e ca talogs of Milky Way 



inere a re zm urivus m tr ie catalogs or Mii Ky way 
HVCs of lPutman et"aLl (pOOl ) and lde Heii eTaTl |2002al) 
with measured fluxes and velocity widths. We assume 
they are distributed with a three-dimensional Gaussian 
distance distribution centered on the Milky Way with a 
given Dhwhm- Given this, we ask for what Dhwhm's 
would we expect zero detections of analogs in our sample 
of six groups. This is done as a Monte Carlo simulation 
with 10,000 trials for a range of Dhwhm between 40-300 
kpc and a population of CHVCs ranging from 27 to 1728 
clouds (0.1-6.4 times the number of cataloged Galactic 
CHVCs). This range in the number of CHVCs per group 
reflects possible variations of the population with group 
mass. If CHVCs are associated with dark matter halos, 
then we expect the total num ber of CHVCs to sc ale in 
proportion to the group mass (jKlypin et al.l ll999) . The 
dynamical ma sses of the se groups are within a dex of the 
Local Group (|Paper HIT ), so this range should be suffi- 
cient. 

We have refined the model for this paper. Previously, a 
CHVC was considered to be detectable if it was above the 
theoretical 10<r sensitivity for its velocity width. Now, 
we use the full completeness function as presented in 
Figure [3] with a cutoff at 5a for each group, acknowl- 
edging that we can detect sources below our previous 
lOer detection limit, but at a significantly decreased level 
of completeness. Effectively, this provides stronger con- 
straints on the distance and mass limits of the CHVCs. 
The cav eats of this model are the same as discussed in 
IPaper J . and fundamentally assume that the properties 
of Milky Way CHVCs are representative of those around 
other galaxies and that our assumed distance distribu- 
tion is a reasonable approximation of the real distribu- 
tion. The former caveat is the most important as our 
survey will not detect the vast majority of CHVCs, but 
only the highest mass ones. This is discussed more below. 

We have also tested a more physically-motivated dis- 
tribution of CHVCs by utilizing a NFW density distri- 
bution: p oc ( r / r3 )(i+( r / ra )ii) > where r s = r vir /C with C 
being the concentration. We vary the value of r s between 
15 and 35 kpc which correspon ds to the rough ran ges of 
C and r v „ for the Milky Way (jKlypin et al.ll2002fh We 



cutoff the dis tribution of CHVCs at lOx r s (150-350 kpc) 
spanning the lKlvpin et all (|002) value of r mr = 258kpc. 
Our choice for this cutoff assumes that all the CHVCs 
are within the virial radius of the Milky Way. 

8.2. Limits on the distances to CHVCs. 

The results of modeling the combined sample of six 
loose groups are shown in Figures 0] and [5] Figure H] 
shows the combined constraints on the number of CHVCs 
per group and Dhwhm of the Gaussian distance distri- 
bution as derived from the non-detection of any CHVC 
analogs in the six groups. Each Dhwhm corresponds 
to an average Mhi for the CHVCs. The figure shows 
that at the 95%, 2a, confidence level a Milky Way- 
like population of 270 CHVCs should be distributed 
with a Dhwhm ^90 kpc implying an average Mhi for 
CHVCs of < 4 x 1O 5 M with the total population hav- 
ing M H i< 1 x 1O 8 M . The median M HI is less than 
1.0 x 1O 5 M and the median CHVC distance is < 116 
kpc. Com pared to our prior limit of Dhwhm ^160 kpc 
(|Paper J ), these limits imply a more tightly clustered, 
less massive population of CHVCs than previously in- 
ferred 5 . The CHVC H I mass function (H I MF) for 
Dhwhm = 90kpc is shown in Figure[6] The shape of the 
H I MF is independent of distance, although the mean 
Mhi will shift with distance. This H I MF is calculated 
statistically for a population of 270 clouds. It is evident 
from this figure that, while we do not expect to detect 
the average CHVC, we expect approximately 10 CHVCs 
to have Mhi> 5 x 10 6 M Q . Our survey is sensitive to 
such clouds, but is highly incomplete for these masses. 
The improvement in our limits as compared to IPaper II 
comes from the refinement of our completeness criteria 
and the addition of three additional groups to our sam- 
ple, although the former has a larger effect than the lat- 
ter. The change in the estimates of the group distances 
had no net effect. 

Figure [5] shows the limits for an NFW halo density 
distribution as a function of r s . At the 95% confidence 
level, with this distribution, the model suggests that r s < 
22 kpc for a Milky Way-like population of 270 CHVCs. 
This corresponds to an average Mhi for CHVCs of < 
3.4 x 10 5 Mq with the total population having Mhi< 
9 x 10 7 M Q . In this model, the median M H i for CHVCs 
is < 4 x 10 4 M Q and the median distance to a CHVC is 
< 90 kpc. These limits are very similar to the Gaussian 
distribution and reflect the uncertainties from assuming 
different distance distributions. It is worth noting that 
a value of r s — 22 kpc corresponds to the best fit model 
param eters for the Milky Way halo from iKlypin et al.l 
(|2002t l with r mr = 258 kpc and C = 12. 

Our model, while simple, appears to be fairly robust 
when considering the possible unique nature of the Milky 
Way CHVC population. If we ignore the most easily de- 
tectable CHVC analogs, our limits are loosened but not 
dramatically. Ignoring the two most easily detectable 
CHVCs in each simulat ion, our distance limits slip to 
match those in IPaper J . The Dhwhm limits are not 
strongly dependent on the total population of CHVCs; 
the limit only decreases by a factor of two when the pop- 



5 In IPaper fl , our derived Mhi limits were incorrect. For 
IPaper ll 's Dhwhm limit of 160 kpc, the average Mhi should be 
< 1O 6 M . 



ulation increases by a factor of 64. As a result, if we were 
to consider all HVCs, not just CHVCs, our limit should 
only decrease to 60-70 kpc. Finally, we can assume that 
the H I fluxes of the 270 CHVCs follow a power law with 
a slope of —2.1 over a range of 4 — 1000 Jy kms™ 1 as 
described in IPutman et~aT1 (|2002f ) and their FWHM ve- 
locity widths are described by a Gaussian with a mean 
of 36 kms" 1 and a disp ersion of 12 kms _1 as described 
bv Ide Heii et alj (j2002al) . The flux limits span the full 
ran ge from the minimum to the maximum cloud flux in 
the IPutman et all (|2002f ) catalog. This removes any ef- 
fects of the discreteness in the cataloged properties on 
the model results. In the end, however, the resulting dis- 
tances and masses are the same with Dhwhm < 100 
kpc. 

For the NFW distribution, there is a stronger depen- 
dence on the total population, but the critical assump- 
tion is the choice of cutoff radius; in this case, we have 
chosen a value similar to the virial radius. If there is 
no cutoff to the distribution, then our model provide no 
distance constraints. Neither model places an artificial 
upper limit on the most massive CHVCs, but the mass 
distribution in both cases truncates around 10 7 M Q . 

The limits are not strongly sensitive to the distance 
estimates to the groups. Dhwhm varies linearly with 
the distance to the group, so that the 20% uncertainty 
in the group distance results in an uncertainty of 20% on 
our derived Dhwhm limits. 

8.3. Comparison of Distance Limits to Past Work 

Our results represent the tightest limits on the distri- 
bution of CHVCs based on observations of galaxy groups, 
but they are roughly consistent wi th the observa tional 
and theoretical results of others. Zwaan (2001) con- 
ducted a similar survey and also failed to detect any 
HVC analogs with resulting mass and distance limits of 
- 10 6 M Q and - 200kpc. Theoretical models o f HVC s 
in dark matter halos bv iMalonev fc Putmanl |2003); 
ISternberg et~ail (|2002f ); lde Heii et all (|2002bf) all suggest 
that the CHVCs should be within ~ 150kpc. Theoret- 
ical models of HVCs as cooling, condensing clouds in a 
hot Galactic halo also place HVCs at distances of ~ 150 
kpc (jMaller fc Bullocki 120041 ) . although detailed simula- 
tions suggest they may be much closer ~ 10 — 60 kpc 
(jKaufmann et alJ l200a ISommer-Larsenl l2006h . Studies 
of HVCs around other galaxies have identified popula- 
tions of clouds which are often associated with the op- 
tical disk, which are probably a s socia t ed with a galac- 
tic fountain (e.g. ISwaters et alj 119971; iBoomsma et al.l 
I2005t iFraternali et alJ 120021 IBarbieri et al.ll2005fh Some 
of these clouds are seen beyond the optical disk, but 
are still within ~ 50 kpc of the galaxy. Other 
H I clouds tha t may have a tidal origin are within 
5 - 100 kpc ([Schneider et all [19831 : ISchneiderl 119851: 



Wilcots et al.lll 996: Hu nter et al.lll99 8: Wilcots & Miller 
1998tlPisano et al.ll 19981 l2002t I Wilcots fc Prescottli2 004) . 
Finally, tho se clouds seen around M 31 are all within 
50 kpc (e.g. iThilker et aD|2004j; |Braun fc Thilkerl 120041: 
IWestmeier et al.l 120051 : iMiller fc Bre gman 200jJ. 

There are also constraints on the distances to Galactic 
HVCs from studies of the population around the Milky 
Way. The only direct constraints come fro m observation s 
of absorption lines towards halo stars (e.g. Wakkcr 2001). 
At present, there are only two complexes with brack- 



eted distances. Complex A is between 8 kpc and 10 kpc 
(|van Woerden et al.Hl99H : IWakker et al.ll2003D. ari d com- 
plex WB is at approximately 8 kpc ( Thorn et al.l l2006). 
However, there are many more distances estimates us- 
ing indirect methods. Many CHVCs show signatures of 
ram-pressure interaction with an ambient medium, such 
as a head-tail or bow shock morphology (|Briins et al.l 
120001: IWestmeier et al.ll2005[) . The required density of 
this medium su ggests that CHVCs should have distances 
~ 50 - 150kpc (iQuilis fc Moordl200~D: iBriins et al.l l200lt 
IWestmeier et al.ll2005D . If HVCs are close enough to the 
Milky Way, then they may be partially ionized by escap- 
ing photons an d emit in Ha. Detec ti ons of Ha emi ssion 
from HVCs bv IPutman et all (f2003h : iWeineri (|200l and 
iTufte et al.l (|2002f) place the clouds within the Galactic 
halo at distances of ~ 40 — lOOkpc. 

8.4. What are the high velocity clouds? 

Our results are independent of the nature of CHVCs. 
They could be associated with a galactic fountain, tidal 
debris, low mass dark matter halos, or condensing clouds 
as long as they are within ~ 90 kpc of the Milky Way, 
have average Mhi of < 4 x 10 5 Af Q , and a total Mhi of < 
10 8 Af Q (considering just the CHVCs). Given these limits, 
can we discriminate between different possible origins for 
CHVCs? 

We have already noted that possible extragalactic 
analogs to HVCs are all observed within ~ 100 kpc of 
the galaxy and completely consistent with our distance 
constraints independent of their proposed origins. As 
seen in Figure El our limits are also consistent with the 
HVCs being distributed like the Milky Way satellites. 
The median distan ce to the Milky Wa y satellites within 
200 kpc is 69 kpc (jGrebel et al.H2003h . Our upper lim- 
its on the median distance to CHVCs are 116 kpc for 
a Gaussian distribution and 90 kpc for an NFW distri- 
bution. If associated with dark matter halos, the Milky 
Way satellites will trace, possibly with some bias, the 
d ark matter distr i bution of the Milky Way. 

iKravtsov et~ail (|2004f ) have identified CDM subhalos 
in their simulations that could accrete gas and form 
stars. The distribution of their simulated luminous satel- 
lites matches that of the Milky Way population of dwarf 
galaxies. The distribution of all dark matter halos has a 
median distance of 116 kpc (jKravtsov et alJ l2004). Fur- 
thermore, they have identified those halos that have asso- 
ciated gas with masses greater than 10 6 M Q ; correspond- 
ing to a Mhi> 10 5 M q . They predict there should be 
50-100 such gas clouds within 200 kpc, and only 2-5 
within 50 kpc. This is fa r fewer than are seen around 
M 31 (jThilker et al.ll2004l ). Given that there are - 300 
CHVCs around the Milky Way, this implies that a sub- 
set of CHVCs could be associated with dark matter halos 
but the majority are not. 

Can the masses of these systems tell us anything 
about the origin of HVCs? For extragalactic H I clouds 
with a putative galactic fountain origin the individual 
clouds have Mjj/~ 10 6 ~ 7 M Q , with the total M H i of all 
such cl ouds around a galaxy ~ 10 7 ~ 9 MrJSchulman et ail 
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2005 



1 19941: iKamphuis fc Sancisil Il993t ISwaters et al 
Fraternali et all 120021 : IBoomsma et all 12004 
IBarbieri et al.l 120051 ) . For tidal debris, they tend to be 
more massive with M H i~ 10 6 ~ 9 M^ dSchneider et al ' 



fl9ll ISchneiderl 119851: IWilcots et all 119961 : iHunter" et al 
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1998t iWilcots fe Milled [19981 : iPisano et all H998I . l200l 
Wilcots fc Prescottl l200l T For HVC analogs associ- 
ated with dark matter halos, mo dels suggest that these 
should have M HI ~ lO 55 " 7 ]^ ([Sternberg etaTI l2002t 
Ide Heii et al"1l2002U iMalonev fc PutmarJl2003t rSimilar 
mass estimates arise from assuming th at HVCs condense 
from the hot gas around galaxies (|Maller fc Bullock] 
[200l . Our limit on the total M m of CHVCs is < 
10 8 M Q . If we were to account f or all of the H VCs, 
then this limit will increase slightly. iPutmanl (2006) sug- 
gest that the total mass (ionized+neutral) in all HVCs 
should be ~ 10 9 M Q if the clouds are within 150 kpc and 
~ 5 x 10 8 M Q if they are within 60 kpc and assuming 
that their neutral fraction is 15%. Our limits imply that 
those galaxies with observed high velocity gas are proba- 
bly not analogous to the Milky Way. In fact most of the 
systems referenced above are either interacting or under- 
going vigorous starbursts. Our limits of the H I mass of 
HVCs can not discriminate between possible formation 
scenarios. 

9. CONCLUSIONS 

We have surveyed six groups of galaxies analogous to 
the Local Group in H I emission using the Parkes tele- 
scope and the ATCA, searching for counterparts to the 
HVCs seen around the Milky Way-H I clouds lacking 
stellar counterparts. No such H I clouds were found. If 
we assume that a Milky Way-like population of HVCs 
is present in each of these groups, we can infer an up- 
per limit on their masses and distances. If HVCs are 
distributed with a three-dimensional Gaussian density 
distribution around galaxies, then the data imply that 
they must be clustered with Dhwhm < 90 kpc, and 
an average Mhi< 4 x 10 5 Af©. While these limits are 
general, we were specifically interested in HVCs poten- 
tially associated with low mass dark matter halos: the 
CHVCs. The total population of CHVCs can only have 
Mhi< 10 8 Af Q with median values D < 116 kpc and 
Mhi< 10 5 M©. This limit is not strongly dependent on 
choosing a Gaussian or NFW distance distribution. Us- 
ing the latter, the median distance drops to 90 kpc with 
the median Mhi< 4 x 10 4 M Q . As such there is not a 
large reservoir of neutral gas around galaxies, and, there- 
fore, any significant reservoir of baryons around galax- 
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TABLE 1 
Group Properties 



Group 


Group Center 8, 
a (J2000) <5 (J2000) 
h:m:s ° : ' : " 


v b 

kms -1 


Distancc c 
Mpc 


Scale 
1° = X kpc 


Num. Galaxies 8 


LGG 93 


03:23:32.1 


-52:13:00 


883 


10.9 


190 


5 


LGG 106 


03:54:56.9 


-47:43:42 


1068 


13.8 


211 


6 


LGG 180 


09:43:54.6 


-31:31:10 


1059 


14.8 


258 


8 


LGG 293 


12:34:12.4 


-07:24:46 


1016 


11.1 


194 


4 


LGG 478 


23:36:12.5 


-36:58:40 


686 


8.6 


150 


4 


HIPASS Group 


13:07:12.4 


-18:31:52 


793 


9.1 


159 


4 



a From Garcia (1993). b Average of individual group galaxies' Vq from NED C Th e distance to the 
group calculated using the multiattractor velocity flow model of Masters (2005) assuming Ho = 
72 km s _1 Mpc -1 . 



TABLE 2 

Parkes Observing Properties 



Group 


Distance 8 
Mpc 


Group Area 
sq. deg. Mpc 2 


Bandwidth 
MHz kms -1 


^0 

MHz 


Beams 


Beamwidth 
14'= X kpc 


mjy 


Sensitivity 15 
10 5 M Q 10 16 


_2 

cm z 


LGG 93 


10.9 


30 


1.1 


8 


1700 


1416 


7 


44 


7.0 


7.2 


3.6 


LGG 106 


13.8 


25 


1.5 


16 


3400 


1414 


13 


56 


5.5 


7.4 


2.S 


LGG 180 


14.8 


25 


1.7 


8 


1700 


1416 


7 


60 


6.5 


11 


3.4 


LGG 293 


11.1 


22 


0.8 


16 


3400 


1414 


13 


45 


6.0 


4.8 


3.1 


LGG 478 


8.6 


35 


0.8 


16 


3400 


1413 


13 


35 


6.5 


3.5 


3.4 


HIPASS Group 


9.1 


48 


1.2 


16 


3400 


1414 


13 


37 


9.0 


5.2 


4.6 


a The distance to the group 


from Table [T] b 


The 1(7 sensitivity over 3.3 kms . 
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Fig. 1. — The number of detected sources as a function of N20 (top) and SNRi„t (bottom) in units of channels and theoretical integrated 
signal-to-noise ratio, respectively. The left two panels show the real sources as filled circles. The right two panels are for the simulated 
sources, with the totals (open squares) and those detected (stars). The open circles in the left panels reflect the real data corrected for 
incompleteness based on the simulated sources in the right panels. No correction is applied if no simulated sources were detected. 
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Fig. 2. — The simulated sources that were detected (solid circles) and not detected (open circles) as a function of integrated flux (in units 
of integrated signal-to-noise ratio) and linewidth (in units of channels) in all six groups. The flux units are defined as the noise in a single 
channel times the square root of the number of channels. For reference, signal-to-noise ratios of 3, 5, and 10 are marked on the figure. 
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Fig. 3. — The completeness for our survey is shown by scaling the HIPASS completeness function from [Zwaari et al.l (2004) to our 
channel size and sensitivity. The abscissa is in units of channels and the ordinate is in units of signal-to-noise in a single channel. The 
solid lines represent (from bottom to top) 50%, 75%, 95%, and 99% completeness. The simulated sources from Figure [2] are plotted here 
for comparison. The dashed lines mark the theoretical 5<r (lower line) and 10cr (upper line) detection limits. 
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<M m > (M Q ) 
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D HWHM ( k P°) 

Fig. 4. — Combined constraints on the Milky Way CHVC population as a function of the number of CHVCs per group and Dhwhm 
(or the average H I mass of the CHVC) assuming a three-dimensional Gaussian distance distribution. The dashed line marks the number 
of CHVCs identified around the Milky Way. 
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Fig. 5.— 
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As in Figure [4] but plotted as a function of r s for an NFW distribution. 
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Fig. 7. — The cumulative distance distribution of Milky Way satellites (solid line with stars), the Kravtsov ct al. (2004) CDM subhalos 
(dotted line), and for a Gaussian distribution with Dhwhm = 90fcpc (solid line) and a NFW distribution with r s = 22kpc (dashed line). 



